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^ ' Abstract 

A new approximate solution of vacuum and stationary Einstein field 
. equations is obtained. This solution is constructed by means of a power 

£SJ ' series expansion of the Ernst potential in terms of two independent and 

dimensionless parameters representing the quadrupole and the angular 
momentum respectively. The main feature of the solution is a suitable 
description of small deviations from spherical symmetry through pertur- 
, bations of the static configuration and the massive multipole structure 

by using those parameters. This quality of the solution might eventually 
provide relevant differences with respect to the description provided by 
^— ^ . the Kerr solution. 
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1 Introduction 



Within the context of axisymmetric Einstein vacuum field equations, the good 
qualities of the Kerr solution pQ are well known and indeed this solution is the 
essential reference in the stationary domain, as is the Schwarzschild solution 
in the static one. We are dealing with an exact solution which is very simple 
and has a very interesting multipolc moment structure in order to describe the 
gravitational field of compact bodies. The Kerr solution is written in terms of 
only one parameter, namely a, which represents the angular momentum (J) per 
unit mass. M being the mass - i.e. the order zero of the multipole moments-, 
all the other moments are proportional to a power of this parameter a equal to 
each multipole order, 0] 

M k = M(ia) k (1) 

It is worthwhile mentioning that all multipole moments of the Kerr solution 
equal the corresponding coefficient appearing in the expansion of its Ernst po- 
tential on the symmetry axis. This characteristic makes it easier to calculate 
the multipole moments of the solution by using the FHP algorithm 0]. 

All these reasons, among others, have led to this solution becoming the one 
most used for a broad range of purposes concerning the physical behaviour of 
test particles or properties of the gravitational fields of celestial sources. 

By means of the multipole moments of any solution one can obtain relevant 
information about the physical properties of the gravitational field, and, for 
instance, the multipole structure Q allows one to estimate the relative impor- 
tance of the successive multipole moments in terms of the parameter a, mainly 
for the sub-extreme case (J << M). 

Nevertheless, since the Kerr solution is defined with only one arbitrary pa- 
rameter, it is not possible to distinguish the different contributions from the 
quadrupole and angular momentum terms, respectively, to any kind of physical 
event described by this space-time. It would be very useful to question how dif- 
ferent a solution is with respect to the Schwarzschild solution in order to know 
which different contributions (or the new physics involved) are introduced by 
a new solution when this has no spherical symmetry and this space-time is no 
longer static. To make such a description, a solution with more than one pa- 
rameter, themselves themselves related to different multipole moments, allows 
one to handle them independently. 

In this sense, some authors [H], |7j have supplied static and stationary solu- 
tions of the Einstein vacuum field equations, with prescribed multipole moments, 
that attempt to describe slight deviations from the space-time with spherical 
symmetry. The Monopole-Quadrupole solution MQ |B| is an exact axisymmetric 
and static solution constructed in power series of a parameter q (dimensionless 
quadrupole moment), whereas the MJ solution 7 is a stationary approximate 
solution describing the space-time of a mass with angular momentum and is 
obtained by an expansion on the symmetry axis of the Ernst potential in power 
series of a dimensionless parameter J . This solution represents successive cor- 
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rections to the spherical symmetry due to the rotation, the sum of the series 
being the pure Monopolc-Dynamic dipole vacuum solution. 

One way to provide the physical content to these solutions consists in estab- 
lishing a link between their multipole moments and quantities measured from 
well defined and physically reasonable experiments. This was the purpose of 
paper |5j , in which we established such a link by calculating the total precession 
per revolution of gyroscopes circumventing the symmetry axis. 

As a second result, the possibility emerges of comparing different axisym- 
metric solutions in term of an observable quantity (angle precessed by the gyro- 
scope). In particular, differences in the contribution of the quadrupole moment 
and the angular momentum for the Erez-Rosen solution |§] and Kerr solution 
£Q were obtained with respect to the MQ^ (see [H], jS] or a brief comment 
below) and MJ solutions respectively. 

Some works ^0] have addressed the description of non-spherical collapse 
and the relevance of this circumstance on the fate of the collapse; and other 
works |2U study the influence of thermal conduction on the evolution of a self- 
graviting system out of hydrostatic equilibrium. Most of the life of a star (at any 
stage of evolution) , may be described on the basis of the quasi-static approxi- 
mation (slowly evolving regime), and this is so, because most relevant processes 
in star interiors take place on time scales that are usually much longer than 
the hydrostatic time scale [25], |2H]. Instead of following the evolution of the 
system a long time after its departure from equilibrium, the system is evalu- 
ated inmediately after such a departure. Here ' immediately' means on a time 
scale of the order of the thermal relaxation time, before the stablishment of the 
steady-state resistive flow. In doing so, it is avoided the introduction of numer- 
ical procedures which might lead to model-dependent conclusions, and for that 
evaluation of the system, an analytic solution becomes fundamental to make use 
of suitable approximations; in particular, it is very important to have a solution 
constructed in a suitable way to describe the contributions to these scenarios 
of both the rotation and the quadrupole deformation separately. On the other 
hand, however, it is only obtained indications about tendency of the object and 
not a complete description of its evolution. Then, it should be clear that, for 
sure, the numerical methods used to solve the complicated mathematical prob- 
lem arising at the study of non-spherical colapse scenarios are much better to 
obtain a complete description. 

As is well knwon, the Weyl metrics ^2] are static axysymmetric solutions to 
vacuum Einstein equations which are given by the line element, 



ds * = - e ^dt 2 + e- M (e 2 ~«(dp 2 + dz 2 ) + p 2 d<\> 2 ) 



(2) 



where metric functions have to satisfy 



(3) 



and 



7, p = P (* 2 p - * 2 2 ) 



7, z = 2p*, p *, z . 



(4) 
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Observe that © is just the Laplace equation for (in the Euclidean space), 
and furthermore it represents the integrability condition for (J2J, implying that 
for any" Newtonian" potential we have a specific Weyl metric. Another inter- 
esting way of writting the general solution of Laplace equation, representing an 
asymptotically flat behaviour, was obtained by Erez-Rosen [§] and Quevedo |12| . 
integrating equations ©, Q in prolate spheroidal coordinates. A sub- family 
of Weyl solutions has been obtained by Gutsunayev and Manko starting 
from the Schwarzschild solution as a seed solution. 

In |23 the MQ solution was presented in the following way 

oo 

*m-q = *,» + + g 2 iv + ••• = £ qaxS i a ' ( 5 ) 

where the zeroth order \I/ g o corresponds to the Schwarzschild solution, and ^ q t 
are series that can be summed to obtain exact solutions of the Gutsunayev- 
Manko family an d the Erez-Rosen family as well ^2] ■ It appears that each 
power in q adds a quadrupole correction to the spherically symmetric solution, 
the first one being referred to as MQW It should be observed that due to 
the linearity of Laplace equation, these corrections give rise to a series of exact 
solutions, and so, the power series of q may be cut at any order, and the partial 
summatory, up to that order, gives an exact solution representing a quadrupolar 
correction to the Schwarzschild solution. 

The existence of so many different (physically distinguishable ) Weyl so- 
lutions gives rise to the question: which among Weyl solutions is better entitled 
to describe small deviations from spherical symmetry?. Although it should be 
obvious that such a question does not have a unique answer ( there is an in- 
finite number of ways of being non-spherical, so to speak), we shall invoke a 
very simple criterion, emerging from Newtonian gravity, in order to choose that 
solution: for example, massive multipole moments of an ellipsoid of rotation, 
with homogeneous density, mass M and axes (a, a, b) read: 

M 2n = K - ' -TVfT-J- , e=(a-b)a, 6 
(2n + l)(2n + 3) 

M 2 „+i = (7) 

because of the factor e™, this equation clearly exhibits the progressive decreasing 
of the relevance of multipole moments as n increases. Same conclusion can be 
obtained for the multipole moments of rotation associated to a homogeneous 
distribution of charge, with total charge Q and ellipsoidal geometry with axes 
(a, a, 6) [Hj: 

■hn = (8) 
(-2)"3Qa 2 "+ 2 e"(l-e/2)"Q 
211+1 ~ (2n+l)(2n + 3)(2n + 5) {) 

being the angular velocity of the rigid rotating distribution. 
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Thus, in order to describe small departures from sphericity, by means of a 
solution of Einstein equations, it would be requiered a solution whose multipolc 
structure shares the property mentioned above. For that reason, the MQ solu- 
tion is particularly suitable for the study of perturbations of spherical symmetry. 
The main argument to support this statement is based on the fact that the pre- 
viously known Weyl metrics (e.g. Gutsunayev-Manko, Manko, 7-metric 
[Tl] . etc) present a drawback when describing quasi-spherical space-times. It 
consists of the fact that its multipole structure is such that all the moments 
higher than the quadrupole are of the same order as the quadrupole moment. 
Instead, as it is intuitively clear and as it is shown above for the classical case, 
the relevance of such multipole moments should decrease as we move from lower 
to higher moments, the quadrupole moment being the moste relevant for a small 
departure from sphericity. 

In the behaviour of geodesies is compared with the spherically symmet- 
ric situation, throwing light on the sensitivity of the trajectories to deviations 
from spherical symmetry. The change of sign in the proper radial acceleration 
of test particles moving radially along the symmetry axis, close to the r = 2M 
surface, and related to the quadrupole moment of the source deserves particular 
attention. 

The ultimate aim of the above works UJ, (see ^jl for general expres- 
sions), and the main contribution of the solution presented here, is to obtain 
a hierarchy of solutions describing pure multipole moment space-times. Since 
1918, when Schwarzschild published its static solution of the Einstein vacuum 
field equations, one of the more active topics of relativistic scientists is devoted 
to obtain exact solutions to those equations. So, up to now a lot of works have 
contributed not only with new solutions but developping several techniques to 
construct new solutions, specifically in the axial symmetry case, nevertheless 
with arbitrary, in general, physical content. Afterwards, main aim aroused was 
to get the correct physical interpretation of the solutions obtained. Some of 
these solutions are already known and studied but in general there is no specific 
method which allows us both to understand which is the physical relevance of a 
new solution as well as to construct solutions with prescribed physical behaviour 
beyond of the symmetry properties of the problem. In [7], a generalization 
of the classical gravitational field to General Relativity theory was planned, 
by describing the solution defined as the sum of different pure multipole con- 
tributions, (following an analogy of the solutions for the classical potential), 
constructed by means of suitable harmonic solutions of the Laplace equation. A 
recent work, [TS] , establishes the reciprocal relationship between static solutions 
with axial symmetry and any prescribed multipole structure. The authors in 
|17| . |T8] show how to obtain the Weyl moments required for constructing the 
pure multipole solutions sought. Although these works are highly relevant, we 
have not yet a generalization of those results to the stationary general case, and 
this is the motivation of the solution presented here. 

The following step to the MQ solution is the pure multipole solution, which 
has only mass, angular momentum and quadrupole moment, referred to here- 
after as the MJQ solution. Here we present an approximation to this solution 
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by means of an expansion of the corresponding Ernst potential. It will be shown 
that the resulting solution is a very good candidate to describe corrections to 
the spherical and static configuration due to the rotation and the quadrupole 
deformation independently. 

For this exterior metric that is presented here (MJQ-solution) it would be 
very interesting to construct a source. The motivation for this is twofold: on the 
one hand, it is always interesting to propose bounded and physically reasonable 
sources of gravitational fields, which may serve as models of compact objects. 
On the other hand, spherical symmetry is a common assumption in the study 
of compact self-gravitating objects (white dwarfs, neutron stars, black holes). 
Therefore it is pertinent to ask, how do small deviations from this assumption 
, related to any kind of perturbation (e.g. fluctuations of the stellar matter, 
external perturbations, etc.), affect the properties of the system?. However, 
for sufficiently strong fields, in order to answer to this question it is necessary 
to deal with non-spherically symmetric solution of Einstein equations. Work 
towards this way is in progress, and, for example, in |21j an interior solution for 
the M QW solution was already obtained. 

The paper is organized as follows: In the next section, all relevant equations 
are given; we develop the construction of the approximate MJQ solution and 
analyze of the behaviour of the solution for limiting cases. Reasons for appre- 
ciating the physical relevance of this solution and its meaningful contribution 
to the description of small deviations from spherical symmetry are given. Fi- 
nally, a discussion of the results is presented in the last section, as well as some 
comments on forthcoming works and the proposals therein. 

2 The MJQ Solution 

2.1 The approximation of the solution 

In 1968 Ernst [221 showed a simplification of the field equations derived from the 
vacuum stationary and axysymmetric Einstein equations. In order to do that, 
Ernst made use of a variational principle on the Lagrange function constructed 
with the metric functions / and lo of the line element describing those metrics, 
which in Weyl canonical coordinates ^2] reads as follows: 

ds 2 = -f(dt - ujdcj)) 2 + f-^e^idp 2 + dz 2 ) + p 2 d<p 2 } (10) 

being /, to and 7 functions only depending on (p, z). 

A very simple way to obtain the Ernst equations, is by means of a suitable 
rewrite of the field equations; a redefinition of the metric function u is made 
by introducing a scalar function W whose existence is guaranteed by one of 
the Einstein field equations 1 . Then, the Ernst potential E is defined as a com- 
plex function whose real and imaginary part are the metric functions / and W 
respectively, the real part representing the norm of the Killing vector. 

J It can be proved that the gradiant of this scalar W is exactly the projection of the twist 
of the time Killing vector, describing stationarity, over the 3-dimensional quotient manifold 
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On symmetry axis, the potential £, which is the transformation £= ^ — — 

from the original Ernst potential, can be expanded by means of a power series 
of the inverse Weyl's canonical coordinate z as follows: 



where p represents the Weyl radial coordinate. 

In |7] a method was proposed for obtaining an approximate solution of the 
stationary vacuum Einstein field equations (MJ-solution or pure monopole- 
dynamic dipole solution) by means of an expansion of the Ernst potential £ in 
power series of a dimensionless parameter, J' , directly related to the angular 
momentum of the solution. The expansion to describe that solution (MJ) is 
quite appropriate, since it is known that its Ernst potential on the symmetry axis 
can be written in terms of successive powers of that parameter. This conclusion 
is obtained from the fact that Fodor, Hoenselaers and Perjes 0] have developed 
an algorithm (FHP) which allows us to calculate the Geroch 3 and Hansen 
2 relativistic multipole moments, related to a vacuum stationary axisymmetric 
solution, in terms of the coefficients (m n ) of the expansion on the symmetry axis 
of the Ernst potential £ (jlip . Both the result obtained up to multipole order 
10 by these authors, as well as the calculations we have carried out up to order 
20 jZj, show that the relation between multipole moments and coefficients m n 
is triangular: that is to say, the multipole moment and the the corresponding 
coefficient m n at every order differ in a certain combination of lesser order rrik 
coefficients. Therefore, these relations can be suitably inverted, and it enable us 
to determine unequivocally the coefficients m n which define the Ernst potential 
£ in terms of the multipole moments for any given solution. For example in 
was obtained the coefficients which characterize the Ernst potential of the 
solution having only massive monopole (the mass) and dynamic dipole (angular 
momentum) . 

We are now looking for the stationary axisymmetric space-time whose mul- 
tipole moment structure only has angular momentum and quadrupole moment 
in addition to the mass. By imposing this condition on the coefficients of the 
series expansion of the Ernst potential along the symmetry axis (lll|) we obtain 
the following expressions (only the fists ones are shown, but they have been 
calculated up to order 20): 



oo 




(11) 



m = M 

mi = iJ 

m 2 = Q 

m 3 = 



m 4 = -MJ 2 + -M 2 Q 

i 8 
m 5 = -—J 3 - i—MJQ 



7 



me - 1m'« + 1m^-|Lm^-§^ ,12) 



J, Q and M being the angular momentum, the quadrupole moment and the 
mass of the solution respectively. This result leads to the Ernst potential £ on 
the symmetry axis of the following type 



oo 

n=0 fc=0 



where J= and q= ^jjp and the double index function F{n,k) is a se- 

ries that contains information about the dependence on the Weyl coordinate, 
z, along the axis with the corresponding numerical factor derived from the 

M 

coefficients m n (I13|) . For instance, it can be seen that F(0, 0) = — and 

M 2 . Z 

.FY 1,0) = — 5-. Furthermore, by considering only k — in the double sum 

1— m z 

of 113f) . then the expression for the MJ solution can be obtained. 

Since we are interested in the description of slight deviations from spherical 
symmetry, we shall take into account the leading terms in the expansion l|13JI 
assuming that q, J << 1, and neglecting terms involving powers or products 
of both parameters. 

Thus, with these considerations the approximate solution up to order q and 
JT, describing the space-time with mass, angular momentum and quadrupole 
moment, is given by the following Ernst potential on the symmetry axis: 



where 





0,z)=£o(p = 0,z)- 






£o(p = 0,z) = 


M 

z 




6(p = o,*) = 


M 2 


F(0,1) 


= 6(p = o,«) = 


OO 



(14) 

(15) 
(16) 



3=1 

with A= — and FS, = — because we know from the FHP 

2J (2j + 3)(2j + l)(2j-l) 

algorithm (see coefficients m n of l|13|l) that the MJQ solution, up to order q 

and J , is defined by the following coefficients: 

m 2j +i = , m 2j = Fg j M 2 S +1 q , Vj > 1 . (18) 

Now, in what follows, and according to the scheme developed in [7|, we 
proceed to solve the Ernst equation |23j . 

(#* - l)Ag = 2C*(V0 a (19) 
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to each order in J and q for the full Ernst potential £, with the following form 

e = & + ixj + (20) 

The set of constants appearing in the general solution of (|19H to each order, will 
be determined by restricting the solution on the symmetry axis and comparison 
with expression (|14p. The zero order in and q gives 

(? 2 - 1)A£ = 2£ (V£ ) 2 , (21) 

and we take £o to be the Ernst potential corresponding to the Schwarzschild 
solution. 

The first order in parameter J , i.e., £1 is already obtained in by solving 
the equation 

- l)A£i - 4£ V£ Va + 2& (V£ ) 2 - , (22) 

and by imposing as a neighbourhood condition the known behaviour of the 
solution on the symmetry axis (TSJ. It should be noted that the approximate 
solution to this order, which in prolate coordinates is written 

t,=io + tiJ =- + ^J , (23) 
x x z 

is the same as that arising from the expansion of the Kerr solution on parameter 

J = i — up to first order. 
M 

The first order in parameter q (£2) must fulfills the following equation: 

(£ 2 - 1)A£ 2 - 4£ V£oV£ 2 + 2£ 2 (V£ ) 2 f- ± | = . (24) 

so 1 

This homogeneous equation can be simplified by means of the following redefi- 
nition of the function £2 : 

A 

sl-l 

which leads to the Laplace equation for the function (2 ■ It can be easily checked 
that, as is known, the equation A^2 = is separable and the general solu- 
tion with regular behaviour on the symmetry axis (y = ±1) as well as in the 
neighbourhood of infinity, affords: 

00 

(2 = J2 h nQn(x)Pn(y) , (26) 
n=0 



C2 = - ( 25 ) 



and therefore: 



„2 00 



^ = ^—^-J2KQ n (x)P n (y), (27) 
n=0 

where {x, y} are prolate coordinates; h%_ are constants of integration, and Q n (x), 
P n {y) are associated Legendre functions of the second kind and Legendre poly- 
nomials respectively. 
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In order to compare the restriction on the symmetry axis of this result (|27[) 
with expression (|17fl . we rewrite (|17l) as follows: 

ll.r2 2 -i 

Up = o,z) = — j— W—b[p = o, z)] , (28) 

z A 2 - 1 

and by carrying out an expansion on the parameter A we have: 

7 2 _ M 2 oo oc 2 oo 

&(p = 0^) = ^EA*5:^]=i/EA^ , (29) 

t=0 j=l " j=l 

3 

with the definition 7? = ■ By using Lemma 4 of the Appendix in [7] we 

»=i 

rewrite this expression in terms of associated Legendre functions of the second 
kind as follows: 

6(P = 0, «) = — II + l)^, 2 ;Q 2i (l/A) , (30) 

72i.2j being the coefficient of the Legendre polynomial of order 2i, with the 
power 2j of the variable. Reordering the sums, we obtain: 

2 _ M 2 00 » 

£ 2 (p = 0,z) = — — ^Q 2i (l/A)(4i + l)^7?i 2i>2j . (31) 

i=i j=i 

We now calculate 7? for our case: 



jq = 5j-(j 2 + 2) = 5 _ 15 1 15 1 
J (2j + 3)(2j + l) 4 8 2j + l 8 2j + 3 



2 



and by comparing expression l|27|) on the symmetry axis (y = ±1, x = — ) with 
H31|) . we conclude that: 

h l n +i =° Vn > (33) 
hg = 

n 

ht n = ~(4n + 1) X] 7 i L 2«.2J Vn > 1 . 
J=i 

Using (|32[1 . it can be verified that in the expression for the constants h\ n , i.e., 
h\ n = -(4n + 1)[V -7 2n 2j - V + y , (34) 

3=0 j=0 J 7=0 J 

for n > 2 the two last terms are equal to zero (see Lemma 2 of the appendix in 
0), and taking into account the orthonormality of Legendre's polynomials we 
have: 

hl n = ~(4n + l) , Vn>2, (35) 
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whereas for n = 1 we see that in (|34[) the last term is not zero (by virtue of 
Lemma 2 in |7j) and therefore: 

1 ^ 

h\ = ~. (36) 

Once we have determined the constants h\ n , we can state that the contri- 
bution on the parameter q to the solution is given by the potential 

C 2 - ~Q*( X ) P *(V) ~ £ ^n+l)Q 2n (x)P 2n (y) . (37) 

71=2 

Alternatively, by using the Heine identity [201 (see also equation (72) in 0) 
we have: 

C 2 = -Qo{x)P (y) - ^Q 2 (x)P 2 (y) - £ -(4n + l)Q 2n (x)P 2n (y) (38) 

= §Qo(x)P (y) - jQ2(x)P 2 (y) ~ ■ 
4 4 41-^ 

2.2 Analysis of the solution 

We shall make some comments on the approximate solution obtained 2 : 



i : = ~ + ^ + i 1 —^liQo^)Po(y) - Q 2 (x)p 2 ( y ) - cw) 



As is well known, the general stationary axisymmetric gravitational field can be 
described by the line element [20] 

ds 2 = -f{dt-ud(t>) 2 + (40) 

a 2 rV 7 (* 2 - y 2 )(^r + 7^) + (* 2 - i)(i - ?/W] , 
— l l — y 

where the coordinate system is {t, x, y, </>}, the space-like coordinates {x, y} be- 
ing the prolate spheroidal coordinates, x represents a radial coordinate, whereas 
the other coordinate y represents the cosine function of the polar angle, with a 
a real constant that, for our case, can be identified with the mass; the unknown 
functions /, 7 and u> only depend on x and y. 

The Ernst potential [201 E = i=| is defined by E = f + iW, where 

W, x = a-^x 2 - 1)- (41) 
W,» = -ff-^l-y 2 )" 1 / 2 ",-. 



2 Note that for a real variable x > 1, we have 1221 the following associated Legendre functions 
of the second kind: Qq{x) = — §lnf=j, Q2<z) = — jl^s) In |=j - fa: 
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and the function 7 satisfies the following partial differential equations: 



1 - v 2 

%* = IW( 2 x ( 42 ) 

Af z {x z - y z ) 



[x{x 2 - l)E tX E* x ~ x(l - y 2 )E, y E* y - y{x 2 - l){E^E* y + E >y E* x ) 



x 2 -I 

7,2/ 



4/2(a ,2_ 2/ 2 ) 

[y(x 2 - l)E jX E% - y(l ~ y 2 )E, y E* y + x(l - y 2 )(E, x E% + E tV E%)] . 
Hence, from l|39f) we obtain: 

/ = ^{l+\q[Q {x)P {y)-Q 2 {x)P 2 {y)-^^}} (43) 
x + 1 2 x — y 

W = ^ - (44) 

(1 + x) 2 M 2 ( } 

and, therefore, by solving equations 1)42(1 and 143(1 up to order q and J we have: 
M x - 1 y ' 



7 = iln(4-4) - ^(1 - y 2 )ln(^-j)+ (-!(,) 
2 x* — y z 8 x + 1 

5 



|?N 2 " Fa 2^(y 2 + ^ - ^ + !)) - 4] ■ 

8 (x z — y^) z 

The solution fulfils asymptotical flatness conditions since all its metric func- 
tions have the good asymptotic behaviour [23j , and the metric does not possess 
any singularity on the symmetry axis, i.e., lim / = 1, lim ui = 0, lim 7 = 0, 

x — >oo x~ >oo x — >oo 

lim 7 = 0. 

It is possible to check the above-mentioned good quality of the solution for 
describing the approximate gravitational held of a mass with only quadrupole 
and angular momentum. First, if we calculate the multipole moments of the 
solution we obtain: 

M = M , Mi = JM 2 = U , M 2 = qM 3 = Q (47) 
M 5 

M 3 = , M 4 = J 2 — 
M 6 

M 5 = —(Jq8-J 3 ) 

M 7 ,„i 51 2 23 

= "21^ IT 9 + TT ^ ' 

These expressions show that the solution possess equatorial symmetry since 
massive multipole moments of odd order are null whereas dynamic multipole 
moments of even order are equal to zero. 
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As can be seen, all multipole moments are several quantities of order higher 
than J, q, or products of both; i.e. the order of approximation considered. 
Take into account that both J and q are very small quantities, since they 
are the dimensionless parameters iJ/M 2 and Q/M 3 respectively, and we are 
considering small deviations from sphericity, i.e., quadrupole moment Q and 
angular momentum J are less than the first multipolar order, the mass M. 
Obviously, this result was expected because the solution is constructed in such 
a way that the coefficients m n are taken up to a linear contribution in both 
parameters. 

Furthermore, this solution has good limits when we consider the parameter 
q or J to be equal to zero. Of course, q = J = reproduces the Schwarzschild 
solution, as we already mentioned, since the solution is constructed by starting 
from £o i the Ernst potential corresponding to spherical symmetry space-time. If 
we take in lp!^|l q = 0, the resulting solution should represent the first correction 
to the static configuration of the Schwarzschild solution due to the effect of 
rotation. In fact, that solution, corresponds to the Kerr solution up to 

order J in the expansion of the Ernst potential in power series of that parameter 
0, i.e., 

^kerr = 1 x^s 

xy/T+J^-Jy 
whose expansion in power series of parameter J gives: 

grr = I (4Q) 

ar 2 x A 

Finally, if we take J = in (|39|l . the corresponding solution 

i = - + q l —^-- A [Q«{x)P»(y) - Q2{x)P 2 {y) - -^-j] (50) 
x x A 4 x z — y z 

should represent the first contribution to the deformation from the spherical 
symmetry configuration due to the quadrupole moment. Therefore, it should 
be the MQ^ solution 0, |Ej constructed with this claim. In order to verify this 
conjecture, it is necessary to check that the metric function ^ corresponding to 
the MQ^ solution generates the Ernst potential (|5U|I up to order q. Since the 
MQ^ solution is given by 'J = + <P&ii we have: 

E=e 2 * =A{l + 2y iq )+(D{q 2 ) , (51) 

with A= e 2 * = s=i, and, 

£ = — = (52) 

l + E B + D v ' 

with the following notation: B = 1 — A, C = 1 + A and D = 2Aq*5>i. Hence, up 
to order q, we have that 

-§-§<§ + U = ; -il^V. + ofcr). (53) 
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By substituting the expression for \ti [5j it is clear that l|53|) is actually equiva- 
lent to the solution obtained (|39ll up to order q, i.e., (|50|l . 

3 Conclusion 

In previous sections an approximate solution of the vacuum Einstein field equa- 
tions with axial symmetry was obtained. The metric functions of the Weyl 
line element are given explicitly by means of the corresponding Ernst potential. 
To construct the solution we have made use of the relationship between the 
coefficients of the Ernst potential appearing at its series expansion in powers 
of the Weyl coordinate along the symmetry axis, and the relativistic multipole 
moments of the solution. 

The relevance of this new solution, in contrast to other known solutions, 
has a very interesting physical meaning since it can be used to describe -in the 
perturbative sense- the corrections to the spherical symmetry that an angular 
and quadrupole moment would incorporate to the solution. As already seen, the 
solution has two independent parameters, q and J , that are directly related to 
these multipole moments and, by using those parameters, the approximation up 
to first order of the pure multipole MJQ solution is constructed. The good limits 
of the solution for the cases q — and J — are verified, thereby recovering 
the MJ solution up to first order (which is exactly Kerr at this order) and the 
MQ' 1 ' solution respectively. 

A significant difference of this solution with respect to the Kerr solution 
is that it allows us to control the magnitude of the contributions from the 
quadrupole deformation and the effect of rotation independently by means of its 
two parameters. The different behaviour of test particles moving along geodesies 
in the radial direction of the MQW solution with respect to that movement in 
the case of spherical collapse is already known JUj . We expect that this solution 
may be a relevant reference to describe non-spherical collapse since it can supply 
more detailed information about the process. Study of all the implications of 
this claim for the solution is still in progress and will be addressed in future 
work, where the good behaviour of the event horizon will be shown as well as a 
detailed study of the test particle geodesies. In any case, it should be stressed 
that the approximation used to construct the solution is correct in the sense 
that it allows one to approach the source without any type of discontinuity, and 
moreso for the strong field case (where the approximation is better since the 
parameters q = Q/M 3 and J = iJ/M 2 become sufficiently less than 1 to be 
neglected). 

Apart from achieving an interior solution for the MQ^ |21|. it is expected 
that a search for an interior solution matching this new solution would provide 
a very interesting global model for the description of compact bodies slightly 
different from the spherical, but never realistic, source. 
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